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Abstract 



Helicity-dependent generalized parton distributions of the nucleon are 
derived from the overlap representation of generalized parton distributions 
using light-cone wave functions obtained in constituent quark models. Re- 
C**) ' suits from two different quark models are used also to study the angular 

momentum sum rule and the spin asymmetry in polarized electron scat- 
tering. 
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1 Introduction 

Over the years hard scattering processes in the deep inelastic scattering (DIS) 
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regime have provided us with considerable insight into the internal nucleon 
structure, i.e. into its quark-gluon substructure as described by Quantum Chro- 
modynamics (QCD). In particular, the spin-independent structure functions 
Fi(x), F2{x), F${x) and the spin-dependent structure functions gi(x), gi{x) 
have been extracted as a function of the fraction x of the quark momentum 
along the direction of the fast moving nucleon (for a collection of data and fit- 
ting programmes, see Ref. pQ). As a consequence, it has been observed that the 
total nucleon momentum and spin are not exhausted by the quark contribution 
alone, and a large debate has originated on the involved quark-gluon dynamics 
(for a recent review, see Ref. [5]). 

A significant step further in the investigation of such a dynamics is offered by 
the recently proposed generalized parton distributions (GPDs) (3J 2] 0D E3 13 E] ■ 
Their importance is due to the fact that they are candidate to give the most 
complete information on the internal nucleon dynamics and to provide us with 
a unifying theoretical background suitable to describe a variety of inclusive and 
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exclusive processes. GPDs are non-diagonal, off-forward hadronic matrix el- 
ements of bilocal products of the light-front quark and gluon field operators. 
As such they carry information on both the longitudinal and transverse dis- 
tribution of partons in a fast moving nucleon and depend on x, the invariant 
momentum square t, and the so-called skewness parameter £ describing the 
longitudinal change of the nucleon momentum. In total, both for quarks and 
gluons there are four distributions conserving the parton helicity and other four 
flipping it 9 . In the case of quarks, in the forward limit GPDs become diag- 
onal matrix elements and three of them reduce to normal quark distributions, 
i.e. the momentum distribution q(x), the helicity distribution Aq(x) and the 
transversity distribution 5q(x). Integrating the quark helicity conserving GPDs 
over x one obtains the nucleon electroweak form factors, that are given in terms 
of off-forward matrix elements of local operators as measured in exclusive reac- 
tions. The second moment of the unpolarized helicity conserving GDPs at t = 
gives a sum rule relating the total quark contribution (including quark orbital 
angular momentum) to the nucleon spin 

GPDs can be probed in deeply virtual Compton scattering (DVCS) and hard 
exclusive production of vector mesons (for recent reviews, see [TU1 ITT1 H2l HI] 
and references therein). First data has become available J3] making the quest 
for modeling GPDs more urgent. 

In the literature there are two approaches used to model the nucleon GPDs. 
One is a phenomcnological construction based on reduction formulae where 
GPDs are related to the usual parton distributions by factorizing the momen- 
tum transfer dependence in agreement with the nucleon electroweak form fac- 
tors p ilfillTTlll7j . This leads to a parameterization of GPDs in terms of double 
distribution functions. However, this approximation is only expected to hold at 
relatively small values of —t I n addition care must be taken to add the 

so-called D-term |19|. an odd function having support < |£| and required to 
satisfy the polynomiality property. This important property follows from Her- 
miticity and time-reversal, parity and Lorentz invariance and implies that the 
m-th moment in a: of a GPD at t — is an even polynomial in £ of degree less 
than or equal to m [TU]. 

Another approach is based on direct calculation of GPDs in specific dynam- 
ical models. The first model calculations were performed using the MIT bag 
model |20j ■ The four helicity conserving GPDs were studied and shown to have 
a quite weak £ dependence, while their t dependence roughly follows the nu- 
cleon form factors behaviour, thus confirming the intuition at the basis of the 
double-distribution assumption. However, the bag model breaks chiral symme- 
try by boundary conditions at the surface, and the initial and final nucleons are 
not good momentum eigenstates. As a consequence, a support violation occurs 
(with GPDs small, but nonvanishing beyond x = 1) and the GPDs behaviour 
in the region \x\ < |£| is not fully reliable. Moreover, an antiquark distribution 
with negative sign is produced when putting three valence quarks in the bag. 

Further calculations have been performed in the chiral quark-soliton model 
The model is based on an effective relativistic quantum field theory 
which was derived from the instanton model of QCD vacuum (see Ref. [22] and 
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references therein). The instanton fluctuations of the gluon field are simulated 
by a pion field binding the constituent quarks inside the nucleon. The model is 
theoretically justified in the limit of the large number of colours N c and satisfies 
all general QCD requirements (sum rules, positivity, inequalities, etc.) including 
the polynomiality property |23|. However, the model comes with an intrinsic 
ultraviolet cutoff in the form of a momentum dependence of the constituent 
quark mass M(p). The effects of the rapidly falling function M(p) are taken into 
account by using some regularization scheme. The uncertainty related to the 
details of the ultraviolet regularization leads to a 10-15% numerical uncertainty 
of the results [23- In the limit of a large number of colours it describes a 
large variety of nucleonic properties typically within 30% accuracy H2| . 
Concerning GPDs, in the leading order in the 1/N C expansion it is not possible 
to obtain results for separate flavours, only special flavour combinations being 
nonzero, i.e. the flavour singlet part of H(x,£,t) and the flavour isovector part 
of the unpolarized E(x, £, t) and of the helicity dependent GPDs, H (x, £, t) and 
E(x,£,t)- 

A complete and exact overlap representation of GPDs has been recently 
worked out within the framework of light-cone quantization [261 12"?] . In a pre- 
liminary investigation of such an approach presented in Ref. |28| a good descrip- 
tion of parton distributions at large x has been achieved with a simple ansatz 
for the wave functions of the three lowest Fock states. In particular, for x > 0.6 
the 80% contribution is produced by the three valence quarks. 

The same approach has been followed recently investigating the link be- 
tween light-cone wave functions (LCWFs) building the overlap representation of 
GPDs and wave functions derived in constituent quark models (CQMs). CQM 
wave functions can be considered as eigenfunctions of the nucleon Hamiltonian 
in the instant-form dynamics and can simply be related to wave functions in 
any form of relativistic Hamiltonian dynamics |30j according to the Bakamjian- 
Thomas construction |31j . Of course, this link is useful in the kinematic range 
where only (valence) quark degrees of freedom are effective. However, in this 
region GPDs are obtained in a covariant approach and exhibit the exact for- 
ward limit reproducing the parton distribution with the correct support and 
automatically fulfilling the particle number and momentum sum rule. 

The method of Ref. was applied to study quark helicity independent 
(unpolarized) GPDs. In this paper we extend it to quark helicity-dependent 
(polarized) GPDs completing the analysis of observable GPDs. In fact, both 
type occurs in DVCS, while hard meson electroproduction is sensitive to un- 
polarized or polarized GPDs depending on whether a longitudinal vector or a 
pseudoscalar meson is produced. In contrast, to date quark helicity flipping 
GPDs seem to contribute only in very peculiar selections of final states in two 
vector meson electroproduction |32) . 

The paper is organized as follows. In Section 2 the relevant definitions are 
summarized and the corresponding expressions in terms of the valence quark 
contribution are given in Section 3. The results obtained within two CQMs 
are presented in Section 4 where they are also applied to study nucleon spin 
asymmetries occurring in inclusive scattering of polarized electrons on polarized 
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targets. Concluding remarks are given in the final section and technical details 
are collected in the Appendix. 

2 The helicity-dependent generalized parton dis- 
tributions 

We work in the so-called symmetric frame |2fill27| . The momentum of the initial 
(final) nucleon is P M (P" 1 ). The average nucleon momentum is then P M = 
I(pM -|_ p'^y The momentum transfer is given by A p = P' M — P M , the invariant 
momentum square is t = A 2 = 2P • A, and the so-called skewness parameter is 
£ = — A + /2P + . We also use the component notation a M = [a + ,aT ,a±] for any 
four-vector with light-cone components = (a°±a 3 )/ y/2 and the transverse 
part o_l = (a 1 , a 2 ). 

The helicity-dependent GPDs are defined starting from the Fourier trans- 
form of the axial vector matrix element 



F q yx (x,U)= ~r~ I dy~ e™ p+ y~ (P> ' ,X'\^y)^ 5 ^y)\P, X) 

47r J v+=v±=o 

(1) 

where the four-vector n is a lightlike vector proportional to (1, 0, 0, —1), A (A') is 
the hclicity of the initial (final) nucleon and the quark-quark correlation function 
is integrated along the light-cone distance y~ at equal light-cone time (y + = 0) 
and at zero transverse separation (y±_ — 0) between the quarks. The resulting 
one-dimensional Fourier integral along the light-cone distance y~ is with respect 
to the quark light-cone momentum k = xP + . The link operator normally 
needed to make the definition (|TJ gauge invariant does not appear because we 
also choose the gauge A + = and assume that one can ignore the recently 
discussed transverse components of the gauge field pHUTI] . 

Following Ref. [S] the leading twist (twist-two) part of this amplitude can 
be parametrized as 



1 

2P" 



F q yx (x,tt) = —u(P',X') 1 + - / 5 u(P,X)H q (x,l; i t) 

■^u(P\X')?^u(P,X)E<i(x,t,t), (2) 
2P 2M 



where u(P, A) is the nucleon Dirac spinor and H q (x, £, t) and E q (x, t) are the 
helicity-dependent GPDs for partons of flavor q, corresponding on the nucleon 
side to an axial-vector and a pseudoscalar transition, respectively. 

An explicit expression of the helicity-dependent GPDs in term of LCWF's 
has been obtained in Refs. |26U27) . Having in mind the link between GPDs and 
CQM wave functions we will restrict our discussion into the region £ < x < 1. 
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In this region and in the symmetric frame 

^ , / -.2-N / N 2-JV 



JV./3 j=l 

x / [c^] N [dk ± ] N 5{x -Xj^y^^r'^x^^eixj)^) 



where Sj labels the quantum numbers of the j-th parton, [3 specifics all other 
quantum numbers necessary for the iV-parton state, and is the helicity of the 
active quark. The set of kinematical variables r,r' are defined as follows: for 
the final struck quark, 

+ ±A+ Xj-£ - ! 1 -Xj -f 

^ = =^T^ = ^ ^=^ + l T3 fA„ (4) 

for the final N — 1 spectators (i ^ j), 

y'i = Yzj^ tf±i = k ±i - ±^-£&±, (5) 

and for the initial struck quark 

^ = =^r7XT = frf = _ T~Tc" (6) 

for the initial N — 1 spectators (i 7^ j), 
Working out the spinor products we have 

F« + {x,S,t) = -F g —(x,t;,t) 

= y/T^H^x,U)--7=^E q (x,^t), (8) 

v 1 -r 



where 

r? = = , (10) 

|A±| 

and 

4£ 2 M 2 , N 
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is the minimal value for — t at given £. 

Using Eq. (J2J, one can derive H q and E q separately from the knowledge of 
Fy X . In particular, E q is directly given by Eq. JJjJ, and 



2M£ 



(12) 



3 The valence-quark contribution 

The valence-quark contribution to GPDs is obtained by specializing Eq. © to 
the case N = 3, i.e. 

F \>x(x,€,t) = — =L=J2J2S S]q sign (fij) J[dx] 3 [dk±\ 3 6(x- Xj) 

» ^ XjTi J=l 

X*#V, {AJ, {t,})*^, {A;}, {r,})6(^), (13) 



where (r, {A^}, {r^}) is the eigenfunction of the light-front Hamiltonian of the 
nucleon, described as a system of three interacting quarks. It is here obtained 
from the corresponding solution \&5^({fcj}, {Aj}, {Tj}) of the eigenvalue equa- 
tion in the instant-form as described in Ref. Separating the spin-isospin 
component from the space part of the wave function, 



*A ] ({*i},{M,M) = ^(*l,fe,te)*Ar(Ai,A 2 ,A3,n,75,73), (14) 



we have 



^\r,{X i },{r i }) = 2(2^ 



1 



-1 1/2 



V>(fcl,fc2,*3) 



_M Zia^s. 

E <A;(^/(fci))^ 2 A ;(i? c/ (fc 2 ))< 2 A i(i? c/ (fc3)) 



X $Ar(A t ljM2:M3,Tl ) r 2 ,T3), 



(15) 



where M is the mass of the non-interacting 3-quark system, = (k^~+k i ) / V2, 
and the Melosh rotations are given by 



Dli 2 (R cf (k)) = (X\R cf (xM ,kM 



. . m + xMq — iff ■ (z x k±) 

(*l ; — Im)- 

(m + xM ) 2 + k\ 



(16) 



In the limit A M — ► 0, where x goes over to the parton momentum fraction 
x, we have 

#«(a?,0,0) = Ag(x), (17) 
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where Aq(x) is the polarized quark distribution of flavor q. Explicitly, the 
following simple expression is obtained 

3 

A.Ti j=l 

x [[^} 3 [dki.} 3 S(x-x j )\^\{x i },{k^i};X i ,T i })\ 2 . (18) 



The polarized quark distribution Aq(x) combined with the unpolarized sin- 
glet quark distribution q(x) — H q (x,0,0) and the unpolarized nonsinglet quark 
distribution E q (x, 0, 0) determines the quark orbital-angular-momentum distri- 
bution L q (x), i.e. PS] 

L q (x) = \{x[q{x) + e*{x)] - Aq(x)}, (19) 

where q(x) = H q (x,0,0) and e q (x) = E q (x,0,0). By integrating Eq. (JTSJ over 
x one recovers the angular momentum sum rule [S] 

J q = i J dxx[q(x) + e q (x)} = iE« + L q , (20) 

where J q is the fraction of the nucleon angular momentum carried by a quark 
of the flavour q, i.e. the sum of spin, 

i£ 9 = J dxAq(x), (21) 

and orbital angular momentum, 

L q = f dxL q (x). (22) 



Thus, from the knowledge of Aq(x) and the unpolarized quark distributions one 
may infer the quark orbital angular momentum L q . 

Furthermore, by integrating the helicity-dependent GPDs over x, for any 
value of the skewness £, one obtains the following relations 

J^dxH q (x,t,t) = G q A (t), J ^dxE q (x,^t)^G q P (t), (23) 

where G q A and G q P are the axial vector form factor and the induced pseudoscalar 
form factor of the quark of flavour q, respectively. 



4 Results and discussion 

In this Section we present results obtained within two constituent quark models, 
i.e. the relativistic version (HYP) [3^1 of the hypercentral quark model [37] 
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Figure 1: The helicity-dependent spin-averaged (H q , left panels) and the helicity-flip 
(E q , right panels) generalized parton distributions calculated in the GBE model for 
flavours u (upper panels) and d (lower panels), at £ = and different values of t: t = 
(solid curves), * = -0.2 (GeV) 2 (dashed curves), * = -0.5 (GeV) 2 (dotted curves). 

and the Goldstone-boson-exchange (GBE) model of Ref. In spite of its 

simplicity the hypercentral model is able to describe the basic features of the 
low-lying nucleon spectrum satisfactorily with a SU(6) symmetric nucleon wave 
function. With the GBE model the baryon spectrum is well reproduced up to 
2 GeV with the correct orderings of the positive- and negative-parity states in 
the light and strange sectors. The resulting nucleon wave functions, without 
any further parameter, yield a remarkably consistent picture of the electroweak 
form factors within a point-form approach to quark dynamics |39| . 

Technical details concerning the derivation of the relevant formulae with 
these models are given in the Appendix. 

The helicity-dependent GPDs calculated in the GBE model for the u and d 
flavours are plotted in Figs. ECU as a function of x at different values of t and 
£. Both H u and H d exhibit a small ^-dependence at £ = 0. At constant t, their 
(rather weak) ^-dependence is entirely governed by the requirement that H q has 
to vanish at the boundaries of the allowed range £ < x < 1 as a consequence 
of including only valence quarks in the present approach. Therefore the peak 
position of H u and H d for increasing £ is shifted to higher values of x. Due 
to the opposite sign of H u and H d in all kinematic conditions, their difference 
is positive and peaked at a value of x comparable to the result obtained in 
the chiral quark-soliton model in the leading order of the 1/N C expansion |18|. 
Also E M and E d have opposite sign as functions of x, their difference being 
rather small at intermediate and large values of a;. In fact, it is known that the 
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Figure 2: The same as in Fig.Elbut for fixed t = -0.2 (GeV) 2 and different values of 
£: £ = (solid curves), £ = 0.1 (dashed curves), £ = 0.2 (dotted curves). 

difference E u — E d is only significantly different from zero in the region \x\ < £, 
where it is dominated by the contribution of the pion pole [18) . 

The helicity-dependent GPDs calculated with the hypercentral and the GBE 
model behave quite similarly (Fig. 0} in spite of the fact that the hypercentral 
model is SU(6) symmetric, while the nucleon GBE wave function contains a 
small SU(6)-breaking part 38, 39 ]. However one should notice that the effect of 
SU(6) symmetry is not equivalent (for spin observables in particular) to the naive 
idea suggested by nonrelativistic dynamics. Within a relativistic approach the 
correlation between motion and spin (helicity) and the large content of high mo- 
mentum components in the wavefunction (determined by the relativistic kinetic 
operator) change the intuitive picture considerably. In fact the SU(6)-breaking 
effects are emphasized (within a relativistic approach) by such correlations and 
high momentum tails, reducing the amount of explicit SU(6)-breaking terms re- 
quired by nonrelativistic approaches. To better appreciate the symmetry effects 
on GPDs, it is worthwhile to discuss first the integral properties of diagonal spin 
and angular momentum observables as defined in Eqs. (20)— (22). 

The different contributions entering the angular momentum sum rule l)2U[> 
are presented in Table Q as obtained in the nonrelativistic SU(6)-symmetric 
quark model, the hypercentral and the GBE model with the same valence quark 
distributions derived here and in our previous paper |29|. The differences be- 
tween the nonrelativistic SU(6)-symmetric model and the hypercentral potential 
are entirely due to relativistic effects. In particular, the unpolarized nonsinglct 
quark distribution e q (x) = E q (x,0,0) vanishes in the nonrelativistic limit when 
all the valence quark are accommodated in the s-wave, while within the rela- 
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Figure 3: The same as in Fig.^but for fixed t = —0.5 (GeV) 2 and different values of 
£: £ = (solid curves), £ = 0.1 (dashed curves), £ = 0.2 (dotted curves). 

tivistic treatment the u and d contributions to e q are quite sizeable in both the 
hypercentral and GBE models. It is remarkable that u quarks contribute much 
more than d quarks. This is expected and in agreement with previous findings 
from lattice QCD (see, e.g., Refs. g0|E]). Only the small difference (0.27 in 
contrast to 0.20) can be ascribed to genuine SU(6)-breaking effects as consid- 
ered by GBE. A similar behaviour can be found for the quark spin contributions 
(E 9 ) and the (quark) angular momentum components (L q ). 

More specifically one can disentangle the effects due to Melosh rotations 
and the high momentum components including (as a "pedagogical" example) 
Melosh rotation into the calculations of spin observables performed with non- 
relativistic wavefunctions (which do not contain high momentum components). 
The effect is somehow surprising: the total quark spin part AE (which reduces 
to 0.46 including u and d components in the HYP model) is enhanced up to 
0.75 021- Let us remark that 0.75 (in contrast to the naive value 1.0 ascribed to 
nonrelativistic quark models) is often quoted as the "the relativistic reduction 
of quark spin contribution to the total nucleon spin due to lower components of 
the wave function" . Our results emphasize that the actual reduction is largely 
influenced by the consistent solution of the mass equation and a simplified guess 
could induce large uncertainties. 

The integral properties shown in the Table 1 share, with all the results of the 
present paper, the limitation of considering valence contribution only. A con- 
sequence is the exact cancellation of the moments J dxxe u (x) + J dxxe d {x). 
Indeed, the constraints of our model, i.e. no gluons, no sea and the momentum 
sum rule exhausted by valence quarks only, lead to J dxx \e u (x) + e d (x)] = 0, 
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Figure 4: The same as in Fig. but for the hypercentral model. 



as shown in Table 1; a result, however, consistent with lattice QCD calcula- 
tions mi- 

In order to exploit all the potential features of the present approach and 
avoid the limitation mentioned, we perform a Leading Order evolution |43| (the 
only available evolution scheme established in the case of angular momentum 
densities) for the moments (20)-(22). One gets jSE] 

iAS(Q 2 ) = iAE(„ 2 ); 

L q {Q 2 ) = (6-50/Bi_i ) i A S( M 2 ) + 6- 50 / 8 %( / x 2 )-^(6- 50 / 81 -l); 

j g (Q 2 ) = &- 50/8 W)-|r 50/81 -i) 

where b = log(Q 2 /A 2 )/ log(^i 2 /A 2 ) and fi 2 = 0.079 GeV 2 and A = 0.232 GeV. 

In the region Q 2 = 1 — 10 GeV 2 , the hypercentral (GBE) model predictions 
range from 0.04 to 0.014 (-0.03 to -0.06) for L q and from 0.20 to 0.25 for J g = 
L g + Ag (and both potential models), results which are in good agreement with 
QCD sum rules predictions ["H] ( J g ~ 0.25) and lattice calculations [IS] (J g = 
0.20 ±0.07). Of course the predictions due to the nonrelativistic models behave 
quite differently, in particular one would expect L q (Q 2 ) ~ — J g ~ —0.25 in the 
range Q 2 = 1 — 10 Gev 2 since AE remains constant in Q 2 . A detailed study [4*2] 
of the quark angular momentum distribution confirms such a behavior showing, 
in particular, a large difference of the results obtained within a relativistic and 
a nonrelativistic approach in the large- a: region. In addition it is important 
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Table 1: Valence contributions to the angular momentum sum rule calcu- 
lated within the nonrelativistic SU(6)-symmetric quark model (NR-SU(6)), the 
SU(6)-symmetric hypercentral model (HYP) and the Goldstone-boson-exchange 
(GBE) model. 





NR-SU(6) 


HYP 


GBE 


J dxxu(x) 


2/3 


2/3 


0.65 


f dxxd(x) 


1/3 


1/3 


0.35 


J dxxe u (x) 





0.20 


0.27 


J dxxe d (x) 





-0.20 


-0.27 


s « 


4/3 


0.61 


0.79 


E' ! 


-1/3 


-0.15 


-0.18 


L" 





0.13 


0.065 


L d 





0.14 


0.13 



to notice that our predictions obtained in the framework of vanishing gluon 
contribution, J ff (/i 2 ) = 0, would not change much considering non vanishing 
Jg(/J?) models. As a matter of fact the scale /i 2 would rise (since at that scale 
the nonvanishing gluons would carry some momentum) and hence b would be 
larger for a given Q 2 . 

Concerning the second moment of GPDs, Eq. I|23(l . we have the following ax- 
ial and pseudoscalar coupling constants: g^ d = G^ d (Q), g P ' d — {M ll /2M)G U p d {t = 
— 0.88M 2 ), with being the muon mass. The calculated values are given in 
Table [3 and are rather far from giving the experimental values. This is evi- 
dently due to neglecting higher Fock states in the present approach and, in the 



Table 2: Axial and pseudoscalar coupling constants calculated with the hyper- 
central (HYP) and the Goldstone-boson-exchange (GBE) model. 





9 u a 


9 d A 


9 u p 


Hp 


GBE 


0.79 


-0.18 


-0.12 


0.04 


HYP 


0.61 


-0.15 


-0.23 


0.07 
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Figure 5: Proton (A\, left panel) and neutron (A™, right panel) spin asymmetries ob- 
tained with the GBE model (upper panels) and the hypercentral model (lower panels) . 
Dashed curves refer to the hadronic scale fig ~ 0.1 GeV 2 , solid curves are the result 
of the evolution at Q 2 = 3 GeV 2 up to NLO. Data from Refs. | 17l H51 H51 ISO] for the 
proton and from Refs. |51l 1521 1331 IF4 55 for the neutron. The predicted SU(6) values 
are indicated by the dot-dashed line. 



pseudoscalar case, to the missing pion-pole contribution dominating the \ow-t 
part of the induced pseudoscalar form factor. 

By studying the relativistic effects introduced by the light-front dynamics 
wc stressed the similarities of the models we are considering, however one could 
also try to investigate better the role of the SU(6)-symmetry breaking intro- 
duced in the GBE model. To this end one can study the spin asymmetry Ai. 
Experimentally, A\ is extracted from the ratio of polarized cross sections in the 
deep inelastic regime as 

A X = ai/2 g3/2 , (24) 

01/2 + °3/2 

where cr 1 / 2 (3/2) is the total virtual photo-absorption cross section for the nucleon 
with a projection of 1/2 (3/2) for the total spin along the direction of the photon 
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momentum. At high Q 2 , A\ can be approximated as the ratio of the spin- 
dependent structure function g\ and the spin-independent structure function 
F<i. In addition, beyond x — 0.3 where mainly valence quarks contribute, A\ 
for protons and neutrons, respectively, is given by 

_ 4Au(x) + Ad(x) Au(x)+4Ad(x) 
Al[X) ~ 4u(x) + d(x) ' Al {X > ~ u{x)+Ad(x) ' (25j 

where u{x) (d(x)) is the unpolarized and Au(x) (Ad(x)) the polarized quark 
distribution of flavour u (d). In the SU(6) limit with both S = and 5=1 
diquark spin states equally contributing, one obtains A\ = 5/9 and A" = 0. 
However, at finite Q 2 one has also to consider the Q 2 dependence of quark 
distributions. This may introduce corrections to the SU(6) limit in spite of the 
fact that the spin asymmetry, being a ratio of distributions, should undergo a 
minimal Q 2 dependence. 

Assuming that the valence quark distributions derived from helicity-indepen- 
dent and helicity-dependent GPDs in the approach discussed here and in our 
previous paper correspond to the hadronic scale /j,q — 0.1 GeV 2 , close to the 
constituent quark mass, Q 2 evolution is introduced along the lines of Ref. |46|. 
Results for A\ and A™ obtained with the GBE and the hypercentral model are 
shown in Fig. [5] In order to compare with available data the hadronic-scale 
results have been evolved to Q 2 = 3 GeV 2 at next-to-leading order (NLO). 
As a consequence of the contribution of sea quarks and gluons introduced by 
evolution, for larger x the results for the proton show an increasing trend shown 
also by data. In the neutron case within CQMs it is known that a relevant 
role is played by SU(6)-symmetry breaking [SSj. Here, the very small SU(6)- 
breaking contribution turns out to be in the opposite direction with respect 
to recent measurements [55] that significantly indicates a positive A\ at x — 
0.6. This discrepancy can ultimately be attributed to the behaviour of the 
unpolarized structure function F% that increases with x (see also Ref. |S7|). Let 
us stress, however, that the approach here described correctly incorporates the 
Pauli principle, a property crucial in understanding SU(6)-breaking effects (see 
Refs. [2302] for discussion). 



5 Conclusions 

The link between light-cone wave functions building the overlap representation 
of generalized parton distributions and wave functions derived in constituent 
quark models was introduced in Ref. (55] along the lines of Refs. 125] 123- R nas 
further been explored here by studying helicity-dependent generalized parton 
distributions of the nucleon. The approach, in principle correct in all kinematic 
regions, has an intrinsic limitation due to the use of the lowest order Fock-space 
components with three valence quark only. Results are therefore obtained in the 
range £ < x < 1 and, when available, can be compared with data only in the 
kinematic range where valence quarks can be assumed as the relevant degrees 
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of freedom. However, they can be useful, e.g., in a Q 2 evolution, and work in 
this direction is in progress. 

Two quark models have been used, i.e. the relativistic hypercentral model of 
Ref. 3BJ and the Goldstone-boson-exchange model of Ref . ■ The results are 
in qualitative agreement with the findings of the chiral quark-soliton model ^H] , 
a modest t dependence is found for H q (x,t;,t) at £ = 0, and the difference 
E u (x,£,t) — E d (x,£,t) is rather small in the valence region. 

In an application to the spin asymmetry in polarized electron scattering the 
same trend of data is obtained for the proton, while the opposite sign at large x 
is found with respect to recent data on the neutron |55| in contrast with previ- 
ous qualitative analyses based on CQMs with SU(6)-symmetry breaking [5B|: a 
difference which can be explained by the correct inclusion of the Pauli principle 
in the present approach. 

As a consequence of the correct relativistic link between wave functions of 
the constituent quark model and light-cone wave functions used in the overlap 
representation of generalized parton distributions, a significant contribution of 
valence quarks to the quark orbital angular momentum is found. This is a 
dynamical relativistic effect leading also to a nonvanishing E(x, £, t) starting 
with only s-wave quarks, a result not obtainable in previous approaches directly 
based on light-cone wave functions. 



A Appendix 

In this appendix we work out the summation over the spin and isospin variables 
appearing in the definition of the helicity amplitudes in Eq. (113(1 . In the case 
of SU(6) symmetric CQM wave functions, the summation over isospin variables 
gives S Tl2 o S T3l/2 + S Tl2 i[S T3 i/ 2 + 2<5 T3 _ 1/2 ]/3 for the proton and c5 Tl2 o <*r 3 -i/2 + 
5ti 2 i[2<5t 3 i/2 + <5 r3 _i/2]/3 for the neutron. On the other hand, the summation 
over the spin variables is carried out in a similar way as in Ref. |29| for the case 
of unpolarized GPDs, by using the explicit expressions of the Melosh-rotation 
matrices appearing in the initial and final light-cone wave function. As a result, 
Eq. (|13fl can be rewritten as 

tl> = 1 ^0^). / n * 4 (» - 1>) s * - *> 

^ i=l \i=l / 

x V 3 {Xyx(* , «) <W + 5*va(« . «)[<W 2 + 2^,-1/2]} . ( 26 ) 



where 



1 LO\LOiLOj, 

M yiy 2 y3 



1p(Hl, K 2 , K3), 



(27) 
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Re 



Im 



Xf + (K,R) = -X°°_(R,k) 

3 

= l[N- 1 (k i )N- 1 {k i )(A 1 A 2 + B x -B 2 )A 3 

3 

X\\{k,k) = -X]}_(k,k) = \{N- l {k l )N- 1 {k l ) 



- (A X A 2 +B 1 -B 2 - iBi, z B2,z)A3 

+2(A 1 B 2rZ + A 2 B ltM )B 3t , 
+2{Bi, x B 2 ^ z + Bi^ z B 2/x )B 3iV 

+2(Bi tV B 2>z + Bi, z B 2 , y )B 3 ^ x , 
Re (xl°_(k',k)' 



Re 



(X°° + (R,R) 



(X°_\(k\k) 



(x^ik.k) 



Im (Xl°_{k ,«) 



{-A X A 2 -BfB 2 + \l > . ., l'.2., )lh. u 

+2{A 1 B 2:X + A 2 B hx )B 3>z 
+2(Bi jX B 2jZ + B\ tZ B 2tX )A 3 
+2(Bi iX B 2<y + B ljV B 2tX )B 3tX , 



lm(xl\(K,K)) = -Im (x^_(k\k)) =Y[N- 1 (K i )N- 1 (K i ) 



i=i 

(-A X A 2 -B!-B 2 + 4B 1>y B 2ty )B 3iX 
+2{AxB 2 , y + A 2 B ltV )B 3tZ 

+ 2{B 1 , X B 2 ,y + Bl,yB 2 ^ X )B 3y y 

+2{B lty B 2>z + B hz B 2 , y )A 3 . 



(28) 



(29) 



]jN-\k[)N- 1 {k i )\{A 1 A 2 +B 1 -B 2 )B 3! y\, (30) 



nTV-i^AT-i^J^^ + Sx-Sa)^], (31) 

i=l 

3 

Re fx} 1 . («,«)) =Y[N- 1 {k[)N- 1 {k l ) 



(32) 



(33) 



In the above equations, N(k), Ai and Bi, with i = 1,2, are defined as in 
Ref. |29| and reported here for convenience 



N{k) = [(m + yMof + k 2 ± } 1/2 . 



(34) 



1G 



Ai = (m + yiMo)(m + yiM ) + Ki iy Ki, y + Ki tX Ki, x , (35) 
Bi, x = -(m + ylM^K^y + im + yiMo)^, (36) 
Bi, y = (m + y' t M^K hX -(m + y t M )K' lx , (37) 

Bi.z — ^i^x^i,y ^i.y^i,x: (^8) 

while A3 and B 3 are given by 

A3 = (m + y 3 M^){m + y 3 M Q ) - n' 3y K 3 ^ y - k 3x k 3iX , (39) 
B 3 , x = (m + y'sM^y + im + yaMoW^y, (40) 
Bs.y = [m + y^M'^K^ + im + ysMoy^, (41) 

B 3,z = K 3 x K 3 , y - K 3; yK 3 , x , (42) 

with yi, yl, and Ki, k^, for i — 1,2,3, defined in Eqs. (@J - 0- 

In the GBE model the nucleon wave functions are expanded on a basis where 
the spin-isospin part is combined with a space part in the form of correlated 
Gaussian functions of the Jacobi coordinates referring to a particular parti- 
tion. The total wave function is a symmetrized linear combination of such basis 
functions over the three possible partitions thus ultimately violating SU(6) sym- 
metry. The calculation with the GBE wave function requires repeating the same 
steps as with the hypercentral wave functions for each partition of the partial 
contribution to the total initial (final) nucleon wave function. 
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